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ABSTRACT 


The quadratic assignment problem (QAP) is a well-known combinatorial 
optimization problem. which could be applied to different applications. The 
main objective of this paper is to present the first discretization of the 
chicken swarm optimization algorithm (CSO) to solve quadratic assignment 
problem without using a local search, the adaptation of CSO in discrete case 
is based on redefining operations and operators of the original version. As 
known, the CSO is a stochastic method inspired from the behavior of 
chickens in swarm while searching for food. The experiments are performed 
on a set of 56 benchmark QAPLIB instances. To prove the choice of the 
adequate parameters, a study is conduct on CSO using simulations on certain 
instances. The discussion of different tests obtains competitive results 


QAPLib compared with the known metaheuristic of Genetic algorithm based on SCX. 
The results demonstrate effectiveness of the proposed CSO-QAP to solve the 
quadratic assignment problem in term of time and quality of solutions. The 
proposed adaptation can be further applied by using a local search strategy 
such as 2-opt in order to solve the same problem or another NP-hard 
combinatorial problem. 


Quadratic assignment problem 
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1. INTRODUCTION 

Combinatorial optimization occupies an important place in discrete mathematics and computer 
science, although the combinatorial optimization problems are often easy to define, but they are generally 
difficult to solve. Indeed, most of these problems belong to the class NP-hard such as Scheduling Workflow 
in Cloud Computing [1], Traveling Salesman Problem [2] and job shop scheduling problem (JSSP) [3]. 
Therefore, NP-hard problems don’t have an effective solution for all the data, then we need to define a formal 
framework for many industry in science, engineering and business. In this work, we intend to discuss one of 
the most interesting combinatorial optimization problem used in the plant layout in order to determine the 
interaction and the distance between two facilities. The aim of the quadratic assignment problem is to fix the 
most effective arrangement of departments within the plant when the flow between departments remains 
constant during the horizontal planning. QAP is a classical NP-hard problem [4] in which it is necessary to 
find the optimal placement of n objects by taking into account both the cost of allocation of an equipment and 
its interaction with other equipments. In the field of location science, many practical problems can be 
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formulated as quadratic assignment problems (QAP) such as the study of Burkard [5] which applied the 
heuristic procedure proposed by Metropolis et al [6] and the approach of Laporte and Mercure [7], which 
used the problem of balancing hydraulic turbine runners. 

Many methods were adapted to solve the quadratic assignment problem, among these adaptations 

are metaheuristics. Metaheuristic is a technique that seems to fit with the structure of any problem. Such as 
simulated annealing by MR Wilhelm in 1987 [8], the taboo search [9], the genetic algorithm [10], the greedy 
genetic algorithm [11], the ant colony algorithm [12] and Hunting Search Algorithm [13]. 
The Swarm optimization algorithms offer a new approach to solve several combinatorial optimization 
problems in engineering and computer sciences. These algorithms are proposed by mimicking the 
inteligence and the behaviors of the population in nature, such as ant colony optimization algonthm 
(ACO) [14], bee colony optimization [15], artificial bee colony algorithm (ABC) [16], bat-inspired 
algorithm [17], Particle Swarm Optimization (PSO) [18], spider monkey [19], firefly algorithms [20] and 
cuckoo search [21]. 

This paper applied the Chicken swarm optimization (CSO) which is an algorithm that simulates the 
behaviors of the chicken swarm. The CSO can effectively harness the intelligence of a chicken swarm to 
solve practically the quadratic assignment problem. The proposed adaptation of CSO 1s the first discretization 
of the chicken swarm optimization algorithm (CSO) to solve quadratic assignment problem, the new 
approach of CSO in discrete case is based on redefinition of operations and operators in order to to solve the 
most of QAPLib instances. 

This paper is organized as follows: In the second section, we describe the quadratic assignment 
problem. In the third section, we present the metaheuristic used in this work to solve the QAP. Then, in the 
fourth section, we provide many simulations to justify the choice of parameters used and the results obtained 
of this contribution. Finally, in the last section, we close this paper with a conclusion. 


2. QUADRATIC ASSIGNMENT PROBLEM 

The management of a plant had long been a field of research, as a result, several approaches are 
developed over the years, but the resource allocation problemis still present. There are different formulations 
to solve the problem of facility management. One of the simple formulations to be modeled is the quadratic 
assignment problem (PAQ). Considered as one of the great challenge in combinatorial optimization, 
Quadratic Assignment Problem is mtroduced in the first time by koopmans and beekmann in 1957 [22]. The 
QAP can be reformulated into many combinatorial optimization problems such as the graph partitioning 
problems and the Travelling Salesman Problem [23]. 

The idea is to set two matrices of size (n,n), given as A = (ajj) which refers to the flows between 


pairs of facilities and B= (Dying »)) for the distance of their locations. After the aforementioned 


descriptions, the purpose of this approach is to find a o* permutation among a set II(n) of permutations, 


which minimizes the cost allocation of facilities mto locations, the mathematical formulation is given as 
follows: 


C(O) = mings Di ae Gybgngc;) (1) 


3. CHICKEN SWARM OPTIMIZATION 

The Chicken swarm optimization (CSO) is a stochastic algorithm inspired by the behaviors of a 
chicken swarm; and introduced by Meng, X.B. And al. [24]. The chicken swarm is divided into several 
groups; each group comprises one dominant rooster and many hens and chicks. Each type of chicken follows 
different movements. 

The roosters, hens and chicks are identified by the fitness value of the chickens themselves. The 
chickens with the best fitness values would be designated to act as roosters or the leader of the group, they 
could search for food in a wider range of space,then they would be followed by the others chickens . 

The chickens with the worst fitness values would be referred as chicks, which could seek for food 
around their mother. The others would be the hens, which randomly choose their groups and establish a 
mother-child relation with chicks. The two relations of dominance and mother-child between the hens and the 
chicks are all updated after few generations ,a parameter G is used to indicate the number of time steps, if 


Ge 2. 20] then the algonthm could give a good results and could easily fall into a local optimum, otherwise 


these statuses remain unchanged. 
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In a combinatorial optimization problem, the position of each chickens represents a solution. The 
roosters with the best fitness values are the most dommants and they could search for food in a wider range 
of space, therefore they have priority to access to food than the hens with lower fitness values, which are 
unstrained to eat food. On the other hand, Chicks move around their mothers to get food. 

In [24] the number of roosters, hens, chicks and mother hens, depicted by RN, HN, CN and MN, 
respectively. The problem is solved in a D-dimensional space, xt;.; stands the positions of all the N chickens 
at time step t while searching for food. 

The position update equation of the rooster can be formulated as: 


xt! =x! #(1+ Randn(0,0°)) 2) 


Lif f<f 
(K-Si) 


( | (3) 
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Where 
- Randn(0,o° )is a Gaussian distribution, O” is a standard deviation. The rooster index k is randomly 


selected from the rooster’s group, and fis the fitness value of the corresponding x. 
The hens can follow their group leaders to forage for food; furthermore, they may randomly locate and steal 
the good food found by other chickens. These situations can be represented as below: 


x = xt 4. S1* Rand * (x, , —x ,)+S2 Rand *(x44,—x,) 


lj 


[ict 
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Where Rand is chosen randomly from [0, 1], rl and r2 are two index/, = V,, the first is the index of the 


rooster otherwise the second is the index of a random chicken from the swarm (rooster or hen). 
At last, the position update equation of the chick is formulated in [24] as follows: 


ae EEX =) (5) 


Where NE [1 N is the index of the chick’s mother and MLE 10, 2 is a randomly selected parameter to 


refer to the relationship between the chicks and its mother. 

As the chicks only learn and get position from their mothers, they will easily fall into the local 
optimum as their mothers. In 2015, Dinghui Wu [25] added to the position update equation of the chicks two 
factors, the first one is the learning factor C which means that the chicks could get information from the 
rooster, the second one is W a Self-learning coefficient for the chicks. The new equation is modified as 
follows: 


PEL 2 t t t t t 
Xj =W*x, | + FL*(X,, j — x )+C#(x1, —x,,) (6) 


Where r is the index of the rooster and mis the index of the mother. 


4. DISCRETE CSO ALGORITHM FOR QUADRATIC ASSIGNMENT PROBLEM 

After observing each individual in the chicken swarm, it was found that the chicken has a capacity 
for communication and ability to learn, thus there is a strict hierarchical order, which can be used to modela 
new stochastic bio-inspired algorithm in order to solve many discrete optimization problem. 

In the begmning, the swarm must be structured by declarng a population, its parameter 
configuration is simple, the fixed identities for each individual is defined after indicating the number of 
roosters, hens and chickens; then the hierarchical order between the identities of the chickens and the 
relations between the mother-child’s are established; And finally, a certanm parameter G must be fixed to 
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regularly update the relations in the chicken swarm so that the problem of falling into a local optimum is 
avoided. 

Each chicken is represented by position, which is surrogated by a vector of N facilities assigned to 
locations, thereby a chicken can search food in a set of solutions S (the search space), and the position of a 
“selected chicken” is the current solution. The purpose of the discretization is to obtain the position of the 
chicken, which provide the minimum of the objective function f:5 >. 

This section describes a new adaptation of the original CSO algorithm; CSO-QAP is established by 
the redefinition of the operators and the operations of CSO (Figure | and 2). 


4.1. Initialize the population 

In this first part, after reading the instances of the QAPLIB, the algorithm begins with the creation 
of an initial population P of N elements randomly selected among p existing permutation of the search space, 
which contains all the possibilities (where p= (n—1)! ). 

2 

Each possibility is obtained by assigning a facility to its location, and then it is represented by an 
array of integer of size n where each integer is an index that refer to the facility’s location. The following 
example is a simple example of a QAP solution that represents an assignment of 6 installations to 6 locations; 
the solution (the table of array) represents a chicken in a known position. 


yp = 
7 1 
tnd 


Location 6 Chicken 2 = (6 {13 |4]5 [2 | = 
Location 4 = J “~~_| After the mouvement | 


Chicken = f4/il3ie6l5]2| 






_-| Before the mouvement 


Chicken1 = al1il3leélsl2” 





Location 3 





Figure 1. Example of Assigning Facilities to Figure 2. Movement of a Chicken in CSO 
Locations in QAP 


4.2. Establish the Hierarchical Order 

The selected population is sorted and ranked after evaluating the fitness value of each individual, 
then a hierarchal order in the swarm is established, the chickens with the minimum fitness values are 
designated as the roosters, the number of the roosters, the hens, the chicks and the mother hens depicted by 
RN, HN, CN and MN. 

The number of roosters can be expressed by the Equation 7, the number of hens by the Equation 8, 
otherwise the number of chicks is calculated by the Equation 9, an index is assigned to each type 1, 2 and 3 
for the roosters, hens and chicks respectively. 


RN = N *rp (7) 
HN =N*hp (8) 
CN = N-(RN+ HAN) (9) 


Where rp is the proportions of roosters and hp is the proportions of hens. MN represents the number of 
mother hens, it could be calculated by 10, and furthermore we assign to each mother a random 1d: 


MN = HN*mp (10) 


Indonesian J Elec Eng & Comp Sci, Vol. 11, No. 3, September 2018 : 925 — 935 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 929 


The sizes of RN, HN, CN and MN directly determine the structure of the swarm, which affect the 
performance of CSO to solve QAP. 


4.3. Create groups 

The swarm is randomly divided into different groups according to the proportions obtained at the 
beginning of the algorithm, the number of groups is equal to number of roosters, we assign rand omly to each 
mother an id hence the relationship between chicks and mother hens is scheduled indiscriminately in 
each group. 

After G iterations, a redistribution of the swarm is required, in order that the roosters find the best 
solution of the swarm and the other individuals could move towards these solutions until the optimum is 
found. 


4.4. The movements of a chicken 

In this adaptation, a movement of a given chicken is represented by the permutation of two of its 
matrix elements. An example of a permutation is described as follows: 

There are three types of chickens (roosters, hens and chicks), each one moves in a structured way 
according to a hierarchy defined in the previous sections. Thus, three types of movements are possible; each 
one is characterized by a set of operations and operators. During the discretization, we defined a set of 
operators, which are @ © and ® as follows: 


4.4.1. Subtraction operation © 
x © x Represents a list of possible permutations applied for solution x; to obtain solution x For 
example : x= 1234567 8andx%=1842765 3 then x © X= (2; 8)(5; 7)(2; 3)(2; 4). 


4.4.2. Multiplication operation © 

This operation means that a random number of permutation is chosen between 0 and | in order to 
apply it to an operation. For example: R=0,5 and x, © x = (2; 8)(5; 7)(2; 3)(2; 4). Then R® (% © x) = 
(2; 8); 7). 


4.4.3. Addition operation © 
This operation indicates that the randomly chosen permutation is applied to move from solution 1 to 
the new solution j. For example: Sl= 12345678 and R® (x © x) =(2; 8) (5; 7) then x= 18347652 


4.4.4. The neighborhood 

In the continuous formula of QAP, the objective is to find a permutation p of n elements that 
minimizes the Equation 1, which uses the fitness value in order to evaluate the solutions. On the other hand, 
in the discrete concept, we can use another approach to say that a solution is close to another neighboring 
permutation, then the differences between the solutions instead of the neighborhood by the fitness value. 


4.4.5. Rooters movement: 

In order to explain why the rooster with the better fitness value can seek for food in a larger space 
,the best roosters will have the opportunity to do more permutations until it find a better solution than the one 
with less fitness value, allowing to search in a wider neighborhood. For that reason, we redefine the operation 
represented in Equation 2 by the Equation 11 as: 


xi" = xi ® Randn(0, 07) ® x! (11) 


In another way, a rooster makes a self-permutation and Randn defines a random percentage of 
permutations between 0 and o to be chosen in order to apply it in the current solution (Rooster ‘position). 


ie : 
Let’s in the example as below x be the current rooster and x is the new solution then 
on 1 
Rand = 20% ,furthermore the current Rooster move to the position 5 a -The value of the parameter Randn 


is proportional to that of oO where Oo = exp(diff) and diff is the percentage of difference between solutions, 
thus the more the differences persist between the current rooster and that of the new position, the more the 
probability of applying permutations increases exponentially, which means that the rooster with better fitness 
values can go to the new position by applying more permutations that is to say search in a wider range of 
space and the value of Randn is chosen randomly between [0,1] ,while for the other roosters of the swarm, 
the value is chosen randomly between 0 and a number less than 1. 
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4.4.6. Hen’s movement 

It is necessary to take mto consideration that hen’s solution must be different to the dominant 
rooster, otherwise in the iterations that follow a sets of roosters will move into the same positions, then the 
search space will be limited. The equation of hens is redefined as follows: 


x =x ®S1®@Randn(x,! x})@S2@Randn®(x,! x) (12) 


The first part of the Equation 12 represents the movement of hens following their leader of group 
with a well-defined probability; the second part means that hens could randomly steal food from other 
chickens (roosters or hens). 


4.4.7. Chicks mouvement 

In this adaptation, the movement of the chicks is divided into 3 steps, firstly the chicks will use a 
self-improvement operation to search for better solutions, then the chicks will look for food by learning from 
their mothers which means that the chicks will have more probability of resemblance between the chicks and 
their mothers. Finally, the chicks will use the same concept to look for the solution of the dominant rooster. 
The equation of hens is redefined as follows: 


XT =WOx OFL@(x! ~)OC@(xi! x) (13) 


Where W is a parameter of self-learning, FL is a learning-Factor, which means that the chick learns from its 
mothers and Cis a learning factor a parameter of resemblance between the chicks and their dominant rooster, 
which means that the chick could learn from of the roosters too. 


4.5. Discrete CSO Algorithm CSO-QAP 
The CSO-QAP in pseudo-code for optimizing the QAP is summarized as follows: 


Discrete Chicken swarm optimization algorithm 


1. Initialize the size of chicken population 

2. Generate a random population of N chickens 

3. Initialize the algorithm’s parameters: N the number of chickens in the 
swarm, Rand, rl is an index of rooster; r2 is an index of chickens, FL, C 
and w. 

4. Evaluate the fitness values at t=O for each chicken (save the global 
best solution) 

5. Rank the chickens and establish a hierarchal order in the swarm 
6. Randomly divide the swarm into different groups 

7. Determine therelationship bet ween the chicks and the mother hens in 


a group. 
8. Finda new solution by updating the position of each rooster, hens and 


chicks using the new equations. 

9. Update the new solution when it is better than the previous one. 
10. Return to step 5 if Gis reached until the maximum number of 
iterations is checked in. 

11. Return results and visualization 


5. EXPERIMENTAL RESULTS AND DISCUSSION 

The performance of the new algorithm is evaluated by the computational experiments were 
performed on the QAP instances extracted from QAPLIB and tested 20 times in 100 iterations for each 
instance. The algorithm was implemented on a DELL in visual studio 2017 using the programing language 
C# and simulated with Intel(R) Core(TM) 17-6500 U CPU 2.5GHZ (4 CPUs) 2.6 GHz and 16.00 GB of 
RAM and Microsoft Windows 10 Professional (64-bit) operating system. 


5.1. Benchmark instances 

In order to evaluate the performance of our new algorithm, we test it on well-known benchmark 
instances of QAPLIB [26], we choose a set of 56 different instances among the 135 instances from the 
QAPLIB, then we compare its results with those of the existing methods in literature, the size of each 
instance is indicated in the instance name. There are many types of instances interpreting and des cribing 
several real and random problems. The first type is the instances obtained from a practical QAP applications 
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Real-life instances, the second is the randomly generated instances with a structure similar to real life 
instances, the third are Grid- based distance matrix instances and finally those in which the distances are 
calculated by the Manhattan distance. 


5.2. Parametric analysis 

In order to choose the optimum parameter values, we execute several tests on chr12a and bur26a. As 
shown in Figure 3, the optimal size of the chicken swarm should be N=500, which provide a better 
compromise between average of BFS and average of time for each run of the instance bur26a. Otherwise 
G=2 guarantees convergence of the algorithm, it may achieve a good results in mmimum average time as 
represented in Figure 7 of the instance chrl2a, this value ensures the robustness of the algorithm and 
guarantees the redistribution of the swarm. 

The number of roosters, hens and chicks could affect the results; we can observe that the percentage 
of roosters must be greater than the percentage of hens as the percentage of chicks to ensure faster 
convergence. In Figure 4 RN should represent 10% of the population with a percentage of 21% for the hens 
as appears in Figure 5, while the rest of the swarm will be chicks respecting that CN=N-RN-HN, then 
CN=69% in Figure 6. We note that the gap between 63 and 72 gives better results for the average of BFS and 
the average of time. The movement of the roosters preserves the exploration; on the other side, the movement 
of hens and chicks performs the exploitation operations. 

Moreover ,as shows in Figure 10, the average of BFS decreases when the self-learning parameter W 
is between 0.4 to 0,6, then we found that for W=0.5, the chicks could search for the solution in a larger range 
of space which avoid the problem of falling into local optimal. Furthermore FL is a randomly chosen number 
between 0,4 and 1, it allows that the chick could learn from the mother which ensure the robustness of the 
algorithm, in Figure 8 we set the parameter FL=0.4. Eventually the chicks could also learn from the leader of 
the group using a learning factor C for better results C=0.7 in Figure 9 which guarantees the optimal balance 
of intensification and diversification for our metaheuristic. 


Table 1. The Parameters values for CSO-QAP 





Parameters Values 
N (population size) 500 
RN (Number of roosters) 10% 
HN (Number of hens) 21% 
CN (Number of chicks) 69% 
G (Number of tours to update 
2 
the algorithm) 
C (Rooster learning factor) 0,7 
FL (Hens learning factor) 0,4 
W (self-learning factor) 0.5 
eae Save(s) RN in chri2a 
5437000 10300 1800 
5436000 10200 1600 
Siren 10100 Aci) 
a Seagal 10000 semis 
> 5434000 9900 1200 
" 5433000 9800 1000 
9700 800 
5432000 9600 600 
5431000 9500 
400 
9400 
a canes 300 500 /OO 900 1100 3300 ee 
. nae ; i | 9200 0 
ne ee 13 5 7 © 21 13 45 17 19° 21. 23 25 27 29 
Figure 3. The average time while varying the Figure 4. The average time and the average BFS while 
population size varying RN 
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Savg(s) HN in chri2a Err (%) Savg(s) CN in chri2a Tavg(s) 
9750 2 9800 1600 

1,8 

‘ 9750 1400 
9700 . 

1,6 

9700 \ 1200 

9650 1A PPP 

1,2 9650 1000 
9600 1 800 

are 9600 

0,8 600 
9550 0,6 occ 

=I 9550 400 
9500 0,4 
95) 9500 

0,2 200 
9450 e O 9450 O 

10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 40 42 44 46 48 50 52 54 56 58 60 62 64 66 68 70 72 74 76 78 80 

Figure 5. The average time and the percentage of Figure 6. The average time and the average of BFS 


ERR while varying HN while varying CN 





9950 5433000 
9900 5432800 
543 ) 
9850 543 260 
5432400 
% 9800 
= “® 5432200 
o uF 
a 9750 go 
o & 5432000 
= 9700 
< 5431800 
9650 
= 5431600 
9600 5431400 
9550 5431200 
0 10 20 30 40 50 60 oO 0,2 0,4 0,6 0,8 1 
Number of iterations Gin chri2a Learning factor FLin bur26a 


Figure 7. The average of BFS while varying the Figure 8. The average of BFS while varying the 


number of iterations G learning-factor FL 
5433200 9700 
5433000 . 
5432800 “i 
5432600 9660 
_ 5432400 964 
8° 5432200 80 
6 5432000 S 9620 
5431800 9600 
5431600 oa 
5431400 nn 
5431200 1560 
Oo 0,2 0,4 0,6 0,8 1 0 0,2 0,4 0,6 08 1 
Learning factor Cin bur26a Self-learning W in chri2a 


Figure 9. The average of BFS while varying the Figure 10. The average of BFS while varying the 
learning-factor C learning-factor W 


Average time of CSO-QAP and SCX 


40 
2c 
5 
3( 
25 
) 
( 
1 i 
10 
x J 7 
*, « * * x 
rd > > xD ~~ ) 
> © >’ © > ~ 


>’ ~s) <>’ <o 2 ~s <2 oo ss a” Oo ao _ _ . _ P 
Y iY V 14 » ay oP a ~ > & &S tai20a.dat tai20b.dat tai25a.dat tai25b.dat tai30a.dat tai4O0b.dat taiSOa.dat tai5Ob.dat 
*) 2 & mc) 2 2 & 2 ) ) > 


2 


@ CSO-OAF scx SOX own CS0-OAP 


Figure 11. Comparison of the percentage ERR of Figure 12. Comparison of the average time of CSO- 
CSO-QAP and SCX QAP and SCX 
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bur26a 
bur26b 
bur26c 
bur26d 
bur26e 
bur26f 
bur26¢ 
bur26h 
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chrl5a 
chr15b 
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chr20c 
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els19 

escl6a 
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escl6e 
escl6g 
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esc32h 
esc64a 
had12 

had14 

had20 

kra30a 
kra30b 
ste36a 
tail2a 
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tail5a 
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tail7a 
tai20a 
tai20b 
tai25a 
tai25b 
tai30a 
tai30b 
tai35a 
tai35b 
tai40a 
tai40b 
tai5Oa 
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tai80b 
tail00a 
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Table 2. Results Obtained By Applying Cso-Qap to Some Qaplib Instances 


Sopt 


5426670 
3817852 
5426795 
3821225 
5386879 
3782044 
10117172 
7098658 
J552 
9896 
7990 
11098 
14142 
3796 
17212548 
68 
292 
160 
16 
28 
26 
130 
Y; 
6 
438 
116 
1652 
2724 
6922 
88900 
91420 
9526 
224416 
39464925 
388214 
51765268 
491812 
703482 
122455319 
1167256 
344355646 
1818146 
637117113 
2422002 
283315445 
3139370 
637250948 
4938796 
458821517 
7205962 
608215054 
1855928 
13499184 
818415043 
21052466 
1185996137 


5.3. Discussion of results 
We performed several tests on a set of different instances. Moreover, the Table 2 summarizes the 
numerical results of the CSO-QAP algorithm for 56 instances of 20 runs, we can notice that the optimal 
solutions is found in 60% of the instances especially with small size. The results in Table 2 summarizes the 
obtained results by applying the CSO-QAP for 56 instances of QAPLIB over 20 independent runs. The first 
column is the name of instance; the Sopt indicated in the column two present the best Known solution of 
each instance in the QAPLIB documentation, the third column present the best found solution by the CSO- 


BFS 


5426670 
3817852 
5426795 
3821239 
5386879 
3782044 
10117172 
7098658 
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8,79 
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9,55 
7,86 
6,14 
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0,26 
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4,31 
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9,19 
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11,87 
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50,45 
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10,02 
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0,015093155 
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0 

0 
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0 
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0 
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1,932 130437 
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0,066 166326 
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0,2 1048 1236 
0,477694806 
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0,482280256 
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0,34032 1432 
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2,25330 1507 
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columns contains the measures use to perform the quality of the solution which are: the best run time T best, 
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the time average Tavg, the average percentage of error Err in 16,the percentage of the Standard Deviation 
PSD in 15 and the last column is the percentage to get the BFS (Succ). The quality of solution is measured 
by: 1) The percentage of the Standard Deviation as: 


psp = 2. x100 (14) 


avg 


Where sp= 





2) The error percentage is calculated by: (15) 


pre = (CxS) 09 


opt 


5.3. Comparison with other metaheuristic 

We compare the average percentage of error to get the best Known solution in our proposed work 
with the results obtained by applying a simple genetic algorithm using sequential constructive crossover for 
the quadratic assignment problem [27]. In this contribution of Genetic algorithm, SCX is applied to solve the 
quadratic assignment problem without using the local search to improve the results, therefore the GA has 
been run on a computer with Intel(R) Core(TM) 17-3770 CPU @3.40GHz and 8.00 GB RAM. Figure 11 
shows that the proposed CSO-QAP is much better in term of the average percentage of error than SCX in the 
GA contribution. 

Furthermore, in order to prove the robustness of the algorithm, Figure 12 presents a comparison 
between CSO-QAP and SCX [27]. The average of time obtained by the CSO-QAP is almost the same of the 
average time obtained by SCX algorithm in most of QAPLIB instances especially in small instance as shows 
in Figure 10; it is obvious that CSO-QAP method is very effective to solve the QAP instances (tai20a, tai20b, 
tai25a, tai25b, tai30a) in a reasonable time. 


6. CONCLUSION 

In this paper, we applied chicken swarm optimization (CSO), as metaheuristic for solving the QAP 
without using a local search. We compare as well the solution found with the best-known solution, which is 
mtroduced in QAPLIB. The results show that the proposed algorithm have effectively demonstrated the 
ability to solve QAP and obtain a promising results in terms of the quality of solutions and the computing 
time .The obtained results are compared with the SCX in the GA contribution. CSO-QAP takes advantage of 
intensification, diversification approaches, the movement of searching for food could be seen as 
mtensification, and the diversification could be represented by the reorganization of the swarm. 

In future research, we can conduct different comparisons between several metaheuristics by using 
different instances of QAP. We can also add low-level heunistics to solve effectively the quadratic 
assignment problem then we will aim to applicate the improved CSO-QAP algorithm in other combinatorial 
optimization problems, especially for problems with high dimensions. 


REFERENCES 

[1] Sheng-Jun Xue and Wu Wu. “Scheduling workflow in cloud computing based on hybrid particle swarm algorithm”’. 
Indonesian Journal of Electrical Engineering and Computer Science, 10(7):1560—1566, 2012. 

[2] Yabo Luo, Min Liu, Zhongxiao Hao, and Dongbo Liu. “An improved nsga-ii algorithm for multi-objective traveling 
salesman problem’. Indonesian Journal of Electrical Engineering and Computer Science, 12(6):4413-4418, 2014. 

[3] Fatima Sayoti, Mohammed Essaid Riffi, and Halima Labani. “Optimization of makespan in job shop scheduling 
problem by golden ball algorithm”. Indonesian Journal of Electrical Engineering and Computer Science, 4(3):542— 
547, 2016. 

[4] Sartaj Sahni and Teofilo Gonzalez. “P-complete approximation problems.” Journal of the ACM (JACM), 23(3):555— 
565, 1976. 

[5] Ramer E Burkard and Franz Rendl. “A thermodynamically motivated simulation procedure for combinatorial 
optimization problems”. European Journal of Operational Research, 17(2):169—174, 1984. 


Indonesian J Elec Eng & Comp Sci, Vol. 11, No. 3, September 2018 : 925 — 935 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 935 


[6] Nicholas Metropolis, Arianna W Rosenbluth, Marshall N Rosenbluth, Augusta H Teller, and Edward Teller. 
“Equation of state calculations by fast computing machines”. The journal of chemical physics, 21(6):1087—1092, 
1953. 

[7] Gilbert Laporte and Hélene Mercure. “Balancing hydraulic turbine runners: A quadratic assignment problem”. 
European Journal of Operational Research, 35(3):378—381, 1988. 

[8] Mickey R Wilhelm and Thomas L Ward. “Solving quadratic assignment problems by ‘simulated annealing’”’. J7E 
transactions, 19(1):107—119, 1987. 

[9] Eric Taillard. “Robust taboo search for the quadratic assignment problem”. Parallel computing, 17(4-5):443-455, 
1991. 

[10] David M Tate and Alice E Smith. “A genetic approach to the quadratic assignment problem”. Computers & 
Operations Research, 22(1):73—83, 1995. 

[11] Ravindra K Ahuja, James B Orlin, and Ashish Tiwari. “A greedy genetic algorithm for the quadratic assignment 
problem”. Computers & Operations Research, 27(10):917—934, 2000. 

[12] Alfredo Reyes Montero and Abraham S’anchez L’opez. “Ant colony optimization for solving the quadratic 
assignment problem’. In Artificial Intelligence (MICAI), 2015 Fourteenth Mexican International Conference on, 
pages 182-187. IEEE, 2015. 

[13] Amine Agharghor and Mohammed Essaid Riffi. “First adaptation of hunting search algorithm for the quadratic 
assignment problem’. In Europe and MENA Cooperation Advances in Information and Communication 
Technologies, pages 263—267. Springer, 2017. 

[14] Marco Dorigo, Vittorio Maniezzo, and Alberto Colorni. “Ant system: optimization by a colony of cooperating 
agents”. [EEE Transactions on Systems, Man, and Cybernetics, Part B (Cybernetics), 26(1):29—41,1996. 

[15] Panta Lu‘ci’c and Du’san Teodorovi'c. “Computing with bees: attacking complex transportation engineering 
problems”. International Journal on Artificial Intelligence Tools, 12(03):375—394, 2003. 

[16] Dervis Karaboga. “An idea based on honey bee swarm for numerical optimization”. Technical report, Technical 
rep ort-tr06, Erciyes University, engineering faculty, computer engineering department, 2005. 

[17] Xin-She Yang. “A new metaheuristic bat-inspired algorithm. Nature inspired cooperative strategies for optimization 
(NICSO 2010)”, pages 65-74, 2010. 

[18] James Kennedy. “Particle swarm optimization. In Encyclopedia of machine learning’, pages 760—766. Springer, 
2011. 

[19] Sandeep Kumar, Rajani Kumari, and Vivek Kumar Sharma. “Fitness based position update in spider monkey 
optimization algorithm”. Procedia Computer Science, 62:442-449, 2015. 

[20] Minmei Huang, Jiyun Yuan, and Jing Xiao. “An adapted firefly algorithm for product development project 
scheduling with fuzzy activity duration”. Mathematical Problems in Engineering, 2015, 2015. 

[21] Amir Hossein Gandomi, Xin-She Yang, and Amir Hossein Alavi. “Cuckoo search algorithm: a metaheuristic 
ap proach to solve structural optimization problems”. Engineering with computers, 29(1):17—35, 2013. 

[22] Tjalling C Koopmans and Martin Beckmann. “Assignment problems and the location of economic activities”. 
Econometrica: journal of the Econometric Society, pages 53—76, 1957. 

[23] Panos M Pardalos, Henry Wolkowicz, et al. “Quadratic Assignment and Related Problems: DIM ACS Workshop”, 
May 20-21, 1993, volume 16.American M athematical Soc., 1994. 

[24] Xianbing Meng, Yu Liu, Xiaozhi Gao, and Hengzhen Zhang. “A new bio-inspired algorithm: chicken swarm 
optimization”. In International Conference in Swarm Intelligence, pages 86—94. Springer, 2014. 

[25] Dinghui Wu, Fei Kong, Wenzhong Gao, Yanxia Shen, and Zhicheng Ji. “Improved chicken swarm optimization”. In 
Cyber Technology in Automation, Control, and Intelligent Systems (CYBER), 2015 IEEE International Conference 
on, pages 681-686. IEEE, 2015. 

[26] Ratner E Burkard, Stefan E Karisch, and Franz Rendl. “Qaplib—a quadratic assignment problem library”. Journal of 
Global optimization, 10(4):391—403, 1997. 

[27] ZH Ahmed. “A simple genetic algorithm using sequential constructive crossover for the quadratic assignment 
problem’. 2014 


Discrete Chicken Swarm Optimization for the Quadratic Assignment Problems (Soukaina C. Bourki Semlali) 


